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SMOOTHING  AND  APPROXIMATION  OF  MULTIVARIATE  FUNCTIONS 

Grant  number  DAAG29- 77-G-0215,  titled  "Smoothing  and  Approximation  of 
Multivariate  Functions,"  was  made  to  The  University  of  Texas  at  Austin  on 
10  September  1977  and  expired  9  September  1979.  Principal  investigator  was 
E.  W.  Cheney,  Professor  of  Mathematics. 

The  research  results  from  this  grant  are  reported  in  8  technical  papers 
which  have  been  or  will  be  published  in  the  open  scientific  literature.  They 
have  also  been  issued  in  Preprint  form  at  the  Center  for  Numerical  Analysis, 
University  of  Texas.  At  the  end  of  this  report,  technical  summaries  of  these 
papers  will  be  presented. 

The  grant  supported  the  principal  investigator  during  four  months  in 
the  summers  of  1978  and  1979.  In  addition,  the  grant  supported  one  Ph.D.  stu¬ 
dent,  John  Respess,  at  a  rate  of  half-time,  during  two  academic  years.  Mr. 
Respess'  dissertation  will  be  completed  in  1980.  The  grant  provided  a  small 
amount  of  supporting  service  through  the  Center  for  Numerical  Analysis.  Several 
machema ticians  from  abroad  who  worked  with  the  principal  investigator  were 
enabled  to  visit  the  University  of  Texas  with  funds  from  this  grant.  Finally, 
some  travel  funds  from  the  grant  enabled  the  principal  investigator  to  accept 
an  invitation  to  lecture  at  the  Third  Biannual  Dundee  Conference  on  Numerical 
Analysis  (Dundee,  Scotland,  June  1979) . 

The  goal  of  the  research  project  was  to  proceed  as  far  as  possible  in 
the  study  of  a  certain  class  of  practical  approximation  problems.  Namely,  the 
problem  is  to  reduce  the  complexity  of  multivariate  functions  by  representing 
them  precisely  or  approximately  by  combinations  of  univariate  functions.  A 
prototype  problem  is  that  of  finding  a  best  approximation  to  a  function  of  two 
variables,  f(x,y),  by  a  sum  g(x)  +  h(y) .  The  prototype  problem  is  thoroughly 
understood  in  the  case  that  the  functions  involved  are  continuous  and  an  approxi¬ 
mation  ^n  the  uniform  sense  is  needed. 

The  general  problem  of  approximatihg  f(x,y)  by 


m 


(1) 


£  a  (x)  y1  +  Z  b  (y)  x1 
i=0  i=0  1 


with  undetermined  continuous  coefficient  functions  a^x)  and  b^(y)  is  very 
poorly  understood,  and  our  efforts  have  been  directed  towards  it. 
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A  typical  practical  computing  problem  which  leads  to  this  type  of 
approximation  is  the  integral  equation 

1 

(2)  u(x)  =  /  f(x,y)u(y)dy 

0 

in  which  u(x)  is  an  unknown  function  and  the  kernel  function  f(x,y)  is  given. 

The  theory  shows  that  u(x)  can  be  determined  easily  if  f(x,y)  is  first  approxi¬ 
mated  accurately  in  the  form  (1) .  The  theory  encompasses,  of  course,  more 
general  equations  than  (2)  and  more  general  forms  than  (1)  . 

Research  Reports  by  E.  W.  Cheney  Supported  by  the  Grant 

1.  "On  Simultaneous  Chebyshev  Approximation"  (with  G.  M.  Phillips,  J.  H.  McCabe), 
C  A  Report  131,  February  1978.  8  pages.  To  appear  in  Journal  of  Approximation 

Theory  26  (1979). 

This  paper  deals  with  the  problem  of  minimizing  an  expression 

max  /  |f  (x,y)-g(x)  |dy 
x 

when  f  is  given  and  g  ranges  over  a  set  of  continuous  functions.  Several  char¬ 
acterizations  of  the  optimal  function  g  are  given.  These  generalize  theorems 
previously  proved  in  the  literature.  This  paper  also  contains  a  bibliography 
of  30  references  to  works  on  "simultaneous  approximation." 


2.  "On  the  Approximation  of  a  Bivariate  Function  by  Sums  of  Univariate  Functions" 
(with  W.  A.  Light),  CNA  Report  140,  August  1978.  27  pages.  To  appear,  Journal 

of  Approximation  Theory. 

This  paper  offers  a  completely  new  proof  of  convergence  of  the  Diliberto- 
Straus  Algorithm  for  uniform  approximation  of  a  continuous  bivariate  function 
by  sums  of  continuous  univariate  functions.  Also,  a  special  analysis  of  the  dis¬ 
crete  case  (which  is  of  importance  in  scaling  problems  of  matrices)  is  given. 

In  particular,  the  functions  which  have  unique  best  approximations  are  completely 
characterized. 
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3.  ''On  the  Algorithm  of  Diliberto  and  Straus  for  Approximating  bivariate  functions 
by  Univariate  Ones"  (withM.  V.  Golitschek) .  Nome r i o a l  Kune  t i ona 1  Ana  lysis 
and  Optimization  1  (1979)  .  341-3b3.  Also,  (.m  Report  l4l,  August  1978. 


This  paper  gives  results  on  the  efficiency  of  the  Di liberto- St raus 
Algorithm  for  obtaining  best  approximations  of  the  form 

f (x.yl  «»  g(x)  t  h(y) 

using  the  Tchebycheff  norm.  For  example,  three  steps  produce  an  approximation 
within  507,  of  the  ultimate  final  accuracy.  Examples  can  be  constructed,  however, 
showing  arbitrarily  slow  linear  convergence.  In  the  discrete  case,  best  approxi¬ 
mations  possessing  certain  extra,  favorable  properties  are  shown  to  exist. 


4.  "The  Approximation  of  bivariate  Functions  by  Sums  of  Univariate  Ones  Using 
the  LpMetric"  (with  W.  A.  Light,  T.  McCabe,  G.  M.  Phillips')  .  CNA  Report  147, 
December  1978,  16  pages.  Submitted  for  publication. 

This  paper  establishes  the  existence  of  best  approximations 

f(x,y>  g(x)  +  h(y) 

in  an  L^-setting.  Thus  the  norm  employed  is  the  L^-norm,  and  the  functions  g 
and  h  can  be  arbitrary  integrable  functions.  It  is  also  proved,  by  an  example, 
that  the  algoritlim  of  Di liberto-Straus  does  not  produce  best  approximations  in 
the  L^-norm. 

5.  "Multivariate  Approximation  with  Tensor-product  Spaces"  (with  W.  A.  Light). 

CNA  Report  153,  November  1979.  25  pages.  Being  readied  for  publication. 

This  report  establishes  the  existence  of  best  approximations  to  a  function 
f(x,y)  by  functions  of  the  form 

n  m 

<*(x,y)  EijWgjfsl  *  Z  b1(a)h1(t) 

1=1  i=l 

when  (1)  the  functions  g^s')  and  h^(t)  are  given  as  essentially  bounded  measurable 
functions;  (2)  the  functions  a^(t)  and  b^(s)  are  allowed  to  range  over  all  inte¬ 
grable  functions;  and  (3)  the  criterion  for  the  approximation  is  the  L^-norm. 
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Thus,  we  want  to  make  this  measure  of  the  error  small: 

//  |f(x,y)  -  (x,y)  |dx  dy 

6.  "Approximation  Theory  in  Tensor-product  Spaces"  (with  John  Respess)  . 

This  is  on-going  work  connected  with  Mr.  Respess'  Ph.D.  research.  A 

representative  result  states  that  each  continuous  function  of  two  variables, 

n 

s  and  t,  possesses  a  best  approximation  of  the  form  y(s)h(t)  +  £,  x  (t)  g  (s) 

i=l  1 

if  h,  g  ,...,gn  are  continuous,  h  has  no  zeros,  (g, . . ., g  )  is  a  Tchebycheff 
set,  and  the  functions  y,  x^,...,x^  are  free  to  range  over  all  possible  continu 
ous  functions. 

7.  "Approximating  Multivariate  Functions  by  Combinations  of  Univariate  Ones," 

Proceedings  of  the  1979  Army  Numerical  Analysis  and  Computers  Conference , 

El  Paso,  Texas,  February,  1979,  pp.  433-440. 

This  is  an  expository  report,  explaining  several  practical  situations 
in  which  a  multivariate  function  must  be  approximated  by  combinations  of  uni¬ 
variate  ones.  Progress  on  the  existence  of  best  approximations  and  on  the 
algorithmic  aspects  of  this  branch  of  numerical  analysis  is  outlined. 

8.  "Best  Approximation  in  Tensor  Product  Spaces,"  CNA  Report  151,  August  1979, 
8  pages.  To  appear,  Proceedings  of  Dundee  Conference  on  Numerical  Analysis 

This  is  an  expository  account  of  the  problem  of  representing  multivar¬ 
iate  functions  by  combinations  of  univariate  ones.  Applications  are  cited  in 
the  fields  of  numerical  treatment  of  partial  differential  equations,  numerical 
solution  of  integral  equations,  and  preconditioning  of  matrices. 
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